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A TWO ZONE SYNCHROTRON MODEL FOR THE 
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ABSTRACT 

The flux and the spectral index in X-ray energy band from the knots of M87 jet as 
observed by Chandra indicate a possible synchrotron origin but cannot be explained 
by simple one zone models with continuous injection of non-thermal electrons. In this 
letter we propose a two-zone model to explain the observed spectra of the knots of 
M87 jet. We consider the synchrotron emission from a region with tangled magnetic 
field where relativistic non-thermal electrons are continuously injected in from an 
associated acceleration region. The acceleration region is assumed to be compact zone 
possibly around a shock front. A power-law distribution of electron is injected into 
the acceleration region and are accelerated to a maximum energy determined by the 
acceleration time scale and the loss processes. With the present model we are able to 
explain the overall broadband features of the knots of M87 jet. Also the present model 
predicts a change in spectral index at ultraviolet energies and future observations at 
these energies can be used to constrain the parameters involved in the model. 
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1 INTRODUCTION 

M87 is a n e arby Riant elliptical galaxy (distance = 16 
Mpc [Xonrvl (| 199 ih ] ) possessing a one-sided jet with pro- 
jected distance « 2 Kpc and bright in radio, optical and 
X-r ay energies. The jet s tructure is ve r y well studied in ra- 
dio fowen et al' (iQSg'); 'Biretta et al." (1995') ; 'Sparks et a_l.| 
lfl996 )). IRJPcrlman ct al. (2001]j^ Sparks ct al. (1996)) and 
optical (Meisenheimer et al.1 (|l996i ): ISparks et al.1 (| 19961 )) 
bands. Prior to Chandra, M87 jet was not very well stud- 
ied at X-ray energies due to the limited angular resolution 
of earlier X-ray telescopes, Einstein and ROSAT. However, 
Chandra due to its better spatial resolution is able to re- 
solve many fainter knots of M87 jet which are observed 
only in radio and optical bands earlier. Moreover, the po- 
sition of these knots in X-ray energie s are nearly coinciden t 
with their radio/op t ical co unterparts (|Perlman et al.l l|200 j ): 
IPerlman fc WilsonI (|2005|)). Also M87 is the o nly radio- 
galaxy (other than CenA( Sreekumar et al.l lll99£))) w hich is 
detected in GeV-TeV energies I Aharonian et al.h 20031 )). Ini- 
tially it was not well understood wh ether the TeV gamma- 
ray e mission region is close to nucleus! Gco rganopoulos et al.l 
ifioOSi )) or from the knot HST- 1 ( Stawarz et al.l (|2006l ) ) . 
However the detection of M87 by HESS confirmed the 



high e nergy emission from the region close to the 
nucleus (| Aharonian et all l|2006l )). Recent VERITA S detec- 
tion of VHE emission from M87 l|Benbowl (|2008l )) again 
suggests the emission may not be from HST-]|j. Consid- 
ering the fact M87 jet is misaligned to the observer, ex- 
plaining this TeV gamma-ray emission required a modi- 
fied model other than the o n e use d to explain blazar emis- 
sion. iNeronov fc AharonianI (|2007l ) explained TeV gamma- 
ray emission due to radiative cooling of electrons accelerated 
by strong rotation induced ele ctric fields in t he vacuum gaps 
in black hole magnetospheres. iLenain et al.l ( |2008, ) proposed 
a multi-blob model with several plasma blobs moving in the 
large opening angle of the jet formation zone. TeV gamma- 
ray emission is explained as Doppler boosted synchrotron 
self Compton radiation by the blobs moving close to the 
line of sight. 

The radio-to-optical emission from the knots of M87 
jet are quite well accepted as synchrotron emission due to 
cooling of re lativistic non-thermal electrons by the magnetic 
field therein (|Perlman et al.l l|200ll )). The fiux and the spec- 
tral indices at X-ray energies indicates a possible continu- 
ation of synchrotron emission of the radio-to-optical spec- 
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tra with the change in the s p ectral index b e yond optical 
energies l|Marshall et al] l|2002h . lStawarz et alj 120051 ) V 



Simple theoretical models viz 
jection model (here af ter CI) 



continu ous 

jKardashevI (Il962h: 



iGinzburg fc Svrovatskiil ll 19681): [ Heavens fc Meisenheimeij 



Meisenheimer et al. I (ll989lV i a nd one tim e injec - 
tion mode l (IJaffe fc Perolal l|l97a ): iKardashevI l|l962h : 



iPacholczvkl (|1970D ) were unable to explain the observed 
X-ray flux and/or the spectral index. The CI model con- 
siders the synchrotron emission from non-thermal electrons 
injected continuously into a region with tangled magnetic 
field. But the X-ray fl ux predicted by this model is more 
than the observed flux l|Perlman fc WilsonI l|2005l )). In one 
time injection model, a burst of el ectrons are injected a t 
t = and allowed to evolve with Jjaffe fc Perolal (Il973l) ) 
or without pitch angle scattering ( KardashevI ( 196^ ) ; 
IPacholczvkl (|l970l )). The model with pitch angle scattering 
under predicts the X-ray flux and the one without pitch 
angle s cattering fail to pre d ict th e observed X-ray spectral 
index (|Perlman fc WilsonI (|2005l )). However CI and one 
time injection models are able to reproduce the broadband 
spectra of the knots of var i ous AGN detected hy C handra 
l|Sahavanathan et al.l (|2003l ): ISambruna et al] (|2002l )). 

In this letter we propose a two zone model to explain 
the non-thermal emission from the knots of M87 jet. We con- 
sider a situation where a power-law distribution of electrons 
are accelerated by an shock. The initial particle spectrum 
which is accelerated by the shock can be an outcome of 
an earlier acceleration process. The particles accelerated by 
the shock cool via synchrotron radiation in a homogeneous 
magnetic field behind it. The present model is similar to 
a situation where particles accelerated by an external shock 
are advected downstream to be accelerated further by an in- 
ternal shock. Two zone models were in fact used by various 
authors to explain the spectr al and the temporal behaviou r 
of blazars llKirk et all (l998l ): [Bhattacharvva et all l|2005h : 
iLi fc Kusunosd (I2OO0I )). In the next section, the formulation 
of the present model and the parameters involved are dis- 
cussed. In §3 we discuss the results of the fitting using the 
present model and compare the present model with the other 
existing models. 



2 THE TWO ZONE MODEL 

We consider the acceleration of a power-law distribution of 
particles (which may be a relic of past acceleration process) 
at a shock front and cool via synchrotron radiation in a 
homogeneous magnetic field. We treat the present scenario 
as a two zones: one around the shock front where the par- 
ticles are accelerated (AR) and downstream region where 
they loose most of their energy through synchrotron process 
(CR). This model is then used to explain the radio-optical- 
X-radio spectra of the knots in M87 jet. We assume the CR 
to be a spherical blob of radius R with tangled magnetic 
field BcR and AR is assumed to be a very thin region with 
magnetic field Bar- Power law distribution of electrons are 
continuously injected into AR and are accelerated by the 
shock characterized by an acceleration timescale tacc. Par- 
ticles are then accelerated at an rate 1/tacc to a maximum 
energy determined by the loss processes. AR is assumed to 



be compact and the emission from CR mainly contributes 
the overall photon spectrum. 

The kinetic equation governing the evolution of elec- 
trons in AR is given by 



d"/ 



dt 
Where 
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iacc 



n(7,t) 



Q{l)dl = i7o7 "^dl for 7min < 7 < 76 



+ Q(7) 
(1) 

(2) 



Here 7 is the Lorentz factor of the electron, teac is the escape 



timescale and Bar ~ t- uafi- 

^ 67rmc 

Equ ation Jl]) can be solve d analytically using Green's 
function (|Atovan fc AharonianI (|l99 9)) and the electron dis- 
tribution for an energy- independent taaa and tesc at time t 
is given by 



n(7,t)=facc7"'"+'' (1 
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Where a = tacc/tesc and the lower hmit of integration Xo is 
given by 



"ymax 



exp{t/ta 



(4) 



"/max = i/{l3ARtacc) > "fb is the maximum Lorentz factor an 
electron can attain in AR. For t » tacc equation Q can be 
approximated to be 7min as the injection term in equation 
Q vanishes for x < ymin- 

The evolution of the electrons in CR is governed by the 
equation 



dt 
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Here the last term is the injection from AR Qajj(7) 
n{'))/teac and for t » tacc 

QARil) ^ goa7-("+^) fl _ ^ ' " ' 



M/]V(7,76) 
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Girmc ^CR' 



and f3cR 

The distribution of electron at time t in CR from equa- 
tion ((5| is given by 



(7) 



Where To = 7/(1 - 7/3cflt). 

From equation (|3|, it can be shown that the injection 
into CR (for a + 1 > p) is a broken power-law with index 
—p for 7 < 7(, and —{a + 1) for 7 > 75. The synchrotron 
losses in CR introduces an additional break jc in the electron 
spectrum depending upon the age of CR (tobs) and the Bcr- 
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Ic = 



PcRtob 



(8) 



The electron spectrum in CR at tabs can then have two 
different spectral shapes depending on the location of 7^ 
with respect to 75. 

(i) 7c > 7i,: The final spectrum will have two breaks with 
indices 



7 



A''(7, tabs) oc < 7' 



-(a+l) 
'{cc + 2) 



< 7 < 7b 

76 < 7 < 7c 

7t: < 7 < 'ymax 



(ii) 7c < 7i): In this case the indices are 



N{'y,tobs) oc < 7 
I 7 



'(p+i) 



< 7 < 7c 
7c < 7 < 76 

7c < 7 < 'ym.ax 



(9) 



(10) 



The synchrotron emissivity e^, from the resultant elec- 
tron spectrum is then calculated by convoluting N{"/,t) with 
single particle emissivity averaged over an isotropic distri- 
bution of pitch angles P(i^, 7) 
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4^/ 
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(11) 



The predicted spectrum by the above model depends on 
nine parameters, which are go, a, jmin, 75, 'Jm.ax, p, Bcr, 
R and tabs- Here a and p are estimated from the radio-to- 
optical and optical-to- X-ray spectral indices, Qo and Bcr are 
constrained using the observed luminosity and equipartition 
magnetic field. For R we assume the physical sizes measured 
m radio l|Hardeel l| 19821 )). Age of the knot toba is chosen to 
introduce a break in the observed spectrum at optical band 
and 7;, is fitted to reproduce the observed X-ray fiux. 7min 
and 7moa: are used as free parameters and are fixed at 5 and 
10^ 
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Figure 1. The spectral fit of the present model with the ob- 
served fluxes in radio(Perlman et al.(2001)), optical(Perlman et 
al.(2001), Waters & Zepf{2005)) and X-ray ((Perlman & Wil- 
son(2005)). The observed spectra of Knot E is fitted with CI 
model (dashed line). Errors in radio and optical arc too small to 
be seen and hence not shown in the figure. 



3 RESULTS AND DISCUSSION 

We applied the above model to explain the knots D, F, A and 
B of M87 jet and the results of the fitting are shown in Fig- 
ure[T]and the parameters used for the fit are given in Table 1. 
The spectrum of knot E can be explained by simple CI model 
and the parameters we quote corresponds to CI model. We 
did not model kno t C due to significan t diffe rences in X-ray- 
optical properties ljPerlman fc Wilson (|2005l )). 

For all the fits shown in Figure [1] 7c < 7;,. However, 
one can fit the spectrum with 7c > 7;, with proper choices 
of the parameters a, 7;, and 7c. This degeneracy arises due 
to the unavailability of UV spectral index since the present 
model predicts the corresponding particle spectral index as 
— {a + 1) or — (p-|-l) (equation ^ and (jlOp ) depending upon 
the above two conditions. Future observations at these pho- 
ton energies may help in validating the present model and 
also will remove this degeneracy. Also, to obtain a precise 
values for p and a, spectral indices at radio-optical and X- 
ray energies should be known accurately. In general tacc and 
tesc can be energy dependent and in such a situation the so- 
lution (equation ([3])) may differ from its form and the index 
beyond 75 may not be the one discussed above. 



Using the present model we estimated the fiux at GeV- 
TeV energies due to inverse Compton scattering and are 
given in Table 1. For the inverse Compton process the tar- 
get photon can be either synchrotron photon(SSC) or cos- 
mic microwave background radiation(IC/CMBR) or both. 
From the observed superluminal velocities of 2 c to 6c and 
the viewing angle of 10° to 20° (jBiretta et al.l (|l999l ')') we 
found the synchrotron photon energy density dominates over 
the Doppler boosted cosmic microwave background radia- 
tion. Hence the flux quoted in Table 1 are due to SSC pro- 
cess and are below the sensitivity of the present gamma-ray 
telescopes M^G/C llMaiumdar fc et"all (|2005l )). HESS and 
the upcoming experiment GLASTQ. Hence discrimination 
of the models based on GeV-Te V observations as discussed 
in iGeorganopoulos et al.l (|2006t ) for 3C273 cannot be done 
for M87 jets. 

A possible scenario of the present model is where AR 
is a region around an internal shock following an exter- 
nal shock. The electron injection into AR can be the one 
which are already accelerated by an external shock and are 

^ http:/ /www-glast. slac.stanford.edu/ 
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advected dow nstream to be accelerated f urther by the in- 
ternal shock (|Tammi fc DempsevI |20o3); iPope fc Melrose! 
l| 19941 )1. An internal shock description for the kno ts in 
AGN j ets has already been d iscussed in literature (iReesI 
l|l978t ): [srhavanathan fc Misral ([2005)). Alternatively, reac- 
celeration of power law electron distribution by turbulence 
at bo undary shea r laye rs can also be another p ossibl e sce- 
nario l|De Yound (|l986l ):[ Stawarz fc OstrowskH (|2003l )). In- 
clusion of these scenario in its exact form into the present 
model will make it more complex and is beyond the scope 
of the present work. 

IPerlman fc WilsonI (j2005l ) proposed a modified CI 
model where the volume within which particle acceleration 
occurs is energy-dependent. This is expressed in terms of a 
filling factor face which is the ratio between the observed 
flux to the flux predicted by the simple CI model . They 
found face declining with increasing distance from the nu- 
cleus suggesting particle acceleration taking place in larger 
fraction of the jet volume in the inner jet than the outer jet. 
The energy dependence of face also indicates that particle 
acceleration regions occupy a smaller fraction of jet volume 
at higher energies. Even though the model is phenomeno- 
logical, it indicates that the process of high energy emission 
from the knots are as complicated as their physical region. 
However the mechanism responsible for the filling factor is 
not explained. 

iStawarz et all (|2006l ) explained the knot HST-1 of M87 
jet as a region when the reconfinement shock reaches the 
jet axis. They considered at the initial stage of M87 jet, 
the particles expand freely decreasing the pressure rapidly 
than the ambient gas pressure. This will develop (in case of 
M87) a reconfinement shock which reaches the jet axis at a 
location which coincide with that of the knot HST-1. They 
postulate this location as the beginning of HST-1 and while 
its outer parts are identified as stationary refiected shock 
formed when the reconfinement shock reaches the jet axis. 
Also they evaluated the ambient radiation field along the jet 
axis and estimated the TeV gamma-ray emission from HST- 
1 initiated by an outburst experienced at the core. Since 
the reconfinement shock requires an initial free expansion, 
the knots downstream HST-1 cannot be explained by this 
model. 

iFleishmanI (|2006l ) explained the fiattening of non- 
thermal spectra in the ultraviolet and X-ray bands ob- 
served from the knots of M87 and 3C273 jets through diffu- 
sive synchrotron radiation(DSR) in random small-scale mag- 
netic fields. Whereas the synchrotron spectrum from regular 
large-scale magnetic field dominates the spectra at low en- 
ergy band. The DSR spectrum at high energy is oc uj~'^ 
where ui is the observed photon frequency and v is the spec- 
tral index of the random magneti c field assumed to be a 
power-law. iHonda fc Hondal (|2007l ') proposed a filamentary 
jet model to explain the observed X-ray spectral index. In 
their model, the jet comprises magnetic filaments of trans- 
verse size A and particles trapped in this filaments are ac- 
celerated by diffusive shock acceleration. The acceleration of 
the electrons bound to a large filament are controlled by the 
radiative losses before escape from the filament. Whereas 
the electrons trapped in smaller filaments escape via ener- 
gization. A critical scale Ac discriminates between the large 
and small scale filaments. They considered a situation where 
the magnetic field is larger for filaments with larger size and 



found the electron energy peaks when trapped in the fila- 
ment of size Ac. The X-ray spectrum is explained by the 
synchrotron radiation of the electrons accelerated in the fil- 
aments of size A > Ac. However, synchrotron radiation from 
large-scale magnetic field itself can reproduce the observed 
X-ray spectrum (present model) involving less number of 
parameters a nd/or assumptions (§2). 

Recently iLiu fc ShenI l|2007f ) proposed a two zone model 
to explain the observed spectra of the knots of M87 jet. In 
their model electrons are accelerated to relativistic energies 
in acceleration region (AR) and loose most of their ener- 
gies in cooling region (CR) through synchrotron process. 
They considered AR and CR are spatially separated and in- 
troduced a break in the particle spectrum injected in CR 
through the advection of particles from AR to CR. This 
along with the cooling break in CR produce a double bro- 
ken power-law with indices —p, — (p+ 1) and — (p + 2) which 
is then used to fit the observed spectra. However the present 
model assumes AR and CR are co-spatial supporting a more 
physical scenario where electrons accelerated by the shock, 
cools in its vicinity. 



4 CONCLUSION 

The observed radio-optical-X-ray spectra from the knots in 
the jets of the FRI radio galaxy M87 are explained within 
the framework of two zone model. We considered a power- 
law electron distribution which are further accelerated in an 
acceleration region and are injected into a cooling region 
where they lose their energy through synchrotron radiation. 
In its simplest form, the model does not consider any spe- 
cific acceleration process but assumes an energy independent 
acceleration timescale. Future observations of M87 knots in 
UV-to-X-ray photon energies will confirm the present model 
and constrain the parameters involved. 

We explored the possibility of the present model to 
reproduce the X-ray fiux of other FRI galaxies (detected 
by Chandra) which are observed to have lower radio lu- 
minosity and relatively smaller jets when compared with 
FRII galaxies. The X-ray emission from FRI jet is quite 
well accepted to be of synchrotron origin whereas for FRII 
and quasars it may be due t o IC/CMBR. However the lat - 
ter is still under debate (see [Harris fc Krawczvnskil (|2006l ) 
for a review about the X-ray emission from extragalactic 
jets). The X-ray emission from the knots and/or the jets of 
the FRI galaxies viz. 3C 6 6B(Hardcastlc ct al. (20 01)), 3C 
346(IWorrafl fc Birkinshawl (l2005h), CenA faardcastle et all 
(|2006l )') and 3C 296 (|Hardcastle et al] (|2005l ') listed in the 
online catalog of extragalactic X-ray jets XJET[f|, which are 
not explained by synchrotron emission from simple one zone 
models, can be reproduced by the present model. 

The author thanks S. Bhattacharyya, N. Bhatt and M. 
Choudhury for the useful discussions and suggestions. The 
author is grateful to referee E. Perlman for useful comments 
and suggestions. This work has made use of the XJET web- 
site. 



http: / /hea-www. harvard.edu/XJET/ 



5 



Table 1. Model Parameters and Knot Properties 



Knot 




a 


P 


76 


BcR 


^obs 


7c 


R 


FlGeV 


F'sOGeV 


F'lTeV 




















(photons / cm^ /s) 


(photons / cm^ /s) 


(photons / cm^ /s) 


D 


7.5 


1.75 


2.35 


1.8 


9.3 


1.6 


5.7 


12 


1.8 X 10-12 


3.7 X lO-i"! 


2.1 X lO-i'' 


E 


4.7 




2.36 




5.9 


2.5 


9.0 


17 


1.7 X 10-13 


4.1 X 10-15 


4.6 X lO-l'^ 


F 


0.6 


2.3 


2.3 


2.1 


4.5 


4.7 


8.2 


29 


1.1 X 10-12 


2.6 X 10-1'' 


1.0 X lO-i'' 


A 


1.0 


2.45 


2.29 


2.0 


4.7 


3.5 


10.1 


55 


2.3 X 10-11 


5.9 X 10-13 


2.1 X 10-15 


B 


0.9 


2.75 


2.3 


1.4 


4.7 


3.9 


9.0 


50 


1.0 X 10-11 


2.3 X 10-13 


3.8 X 10-1® 



Columns:- 1: Knot name. 2: Normalisation (in units of 10— l^) of the power-law injected into AR(for knot E it is the 
normalisation of the power-law injected into CR (see text)). 3: Ratio between the acceleration timescale and escape 
timescale in AR. 4: Index of the power-law spectrum injected into AR. 5: Maximum energy of the electron Lorentz factor 
injected into AR (in units of 10®). 6: CR magnetic field (in units of IQ-* G). 7: Time of observation (in units of 10® 
sec). 8: Break Lorentz factor of the electrons due to synchrotron cooling in CR (in units of 10""). 9: Size of CR (in parsec) 
measured in radio (Hardee (1982)). 10: Flux at IGeV. 11: Flux at 50GeV. 12: Flux at ITeV. 

For all cases, the minimum Lorentz factor(7mi„)injected into AR is 5 and maximum Lorentz factor(7mai) attained in AR 
is 10*. 
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